Mon. Not. R. Astron. Soc. OOP, [T|fT5l (2011) Printed 20 January 2013 (MN WT^ style file v2.2) 



The effect of a massive object on an expanding universe 



Roshina Nandra ' , Anthony N. Lasenby^' and Michael P. Hobson 

^Astrophysics Group, Cavendish Laboratory, JJ Thomson Avenue, Cambridge CB3 OHE, U.K. 

^ Kavli Institute for Cosmology, c/o Institute of Astronomy, Madingley Road, Cambridge CBS OHA, U.K. 



Accepted — . Received — ; in original form 20 January 2013 



ABSTRACT 

A tetrad-based procedure is presented for solving Einstein's field equations for 
spherically-symmetric systems; this approach was first discussed by Lasenby, Doran 
& Gull in the language of geometric algebra. The method is used to derive metrics 
describing a point mass in a spatially-flat, open and closed expanding universe re- 
spectively. In the spatially-flat case, a simple coordinate transformation relates the 
metric to the corresponding one derived by McVittie. Nonetheless, our use of non- 
comoving ('physical') coordinates greatly facilitates physical interpretation. For the 
open and closed universes, our metrics describe different spacetimes to the corre- 
sponding McVittie metrics and we believe the latter to be incorrect. In the closed 
case, our metric possesses an image mass at the antipodal point of the universe. We 
calculate the geodesic equations for the spatially-flat metric and interpret them. For 
radial motion in the Newtonian limit, the force acting on a test particle consists of the 
usual inwards component due to the central mass and a cosmological component 
proportional to r that is directed outwards (inwards) when the expansion of the uni- 
verse is accelerating (decelerating). For the standard ACDM concordance cosmology, 
the cosmological force reverses direction at about z « 0.67. We also derive an invariant 
fully general- relativistic expression, valid for arbitrary spherically-symmetric systems, 
for the force required to hold a test particle at rest relative to the central point mass. 
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1 INTRODUCTION 

Among the known exact solutions of Einstein's field equa- 
tions in general relativity there are two commonly studied 
metrics that describe spacetime in very different regimes. 
First, the Friedmann-Robertson- Walker (FRW) metric de- 
scribes the expansion of a homogeneous and isotropic uni- 
verse in terms of the scale factor R{t). The FRW metric 
makes no reference to any particular mass points in the 
universe but, rather, describes a continuous, homogeneous 
and isotropic fluid on cosmological scales. Instead of using a 
'physical' (non-comoving) radial coordinate r, it is usually 
written in terms of a comoving radial coordinate f, where 
r = fR(t), such that the spatial coordinates of points mov- 
ing with the Hubble flow do not depend on the cosmic time 
t. Here the comoving coordinate f is dimensionless, whereas 
the scale factor R{t) has units of length. For a spatially-flat 
FRW universe, for example, using physical coordinates, the 
metric is 

ds^ = [1 - r^if^(t)] dt^ + 2rH{t) drdt - dr^ - r^dQ.'^, (1) 
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which becomes 

ds'^ = dt^ -R'^{t){df^ +r^dQ.'^) (2) 

in comoving coordinates, where dQ,^ = d6^ + siv? 6d(jP , 
H{t) = R'{t)/R{t) is the Hubble parameter and primes de- 
note differentiation with respect to the cosmic time t (we 
will adopt natural units thoroughout, so that c = G = 1). 

Second, the Schwarzschild metric describes the spher- 
ically symmetric static gravitational field outside a non- 
rotating spherical mass and can be used to model space- 
time outside a star, planet or black hole. Normally the 
Schwarzschild metric is given in 'physical' coordinates and 
reads 

(3) 

but another common representation of this spacetime uses 
the isotropic radial coordinate r, where r — {l + ^) r, such 
that 

— , o 

l-h^V(df'+r'dn'). (4) 
2,-./ ~ 2r/ 

The main problem with the Schwarzschild metric in a cos- 
mological context is that it ignores the dynamical expanding 
background in which the mass resides. 
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McVittie ( 1933 1956 1 combined the Schwarzschild and 
FRW metrics to produce a new spherically-symmetric met- 
ric that describes a point mass embedded in an expanding 
spatially- flat universe. McVittie demanded that: 

(i) at large distances from the mass the metric is given 
approximately by the FRW metric (|2|; 

(ii) when expansion is ignored, so that R{t) = _Ro, one 
obtains the Schwarzschild metric in isotropic coordinates Q 
(whereby f = ri?o, with r defined below); 

(iii) the metric is a consistent solution to Einstein's field 
equations with a perfect fluid energy-momentum tensor; 

(iv) there is no radial matter infall. 

McVittie derived a metric satisfying these criteria for a 
spatially-fiat background universe; 



ds 



2rR(t) 



2rR.{t) 



m \ ' 
2rR{t) ) 



R^{t){dr'^ +r'^dn^), (5) 



where r has been used to indicate McVittie's dimensionless 
radial coordinate, rather than our 'physical' coordinate. In 
Section |3.1| we will see how the two are related and point 
out some of the problems with (fsl. One sees that Q is a 
natural combination of Q and Q; nonetheless, there has 
been a long debate about its physical interpretation. This 
uncertainty has recently been resolved by Kalopcr, Kleban| 
& Martini (120101) andlLake & Abdelqaderl (120111), who have 



shown that McVittie's metric does indeed describe a point- 
mass in an otherwise spatially-fiat FRW universe. In this 
paper we have also independently arrived at the same con- 
clusion, as discussed in Section [3?T| McVittie also generalised 
his solution to accommodate spatially-curved cosmologies, 
which are discussed further in Section [S] 

For a given matter energy-momentum tensor, Einstein's 
field equations for the metric g^j^ constitute a set of non- 
linear differential equations that are notoriously difficult to 
solve. Moreover, the freedom to use different coordinate sys- 
tems, as illustrated above, can obscure the interpretation of 
the physical quantities. In a previous paper, |Lasenby et al.| 
( 1998 1 presented a new approach to solving the field equa- 



tions. In this method, one begins by postulating a tetrad (or 
frame) field consistent with spherical symmetry but with un- 
known coefficients, and the field equations are instead solved 
for these coefficients. In this paper, we follow the approach 



of Lasenby et al. ( 1998 1 to derive afresh the metric for a 
point mass embedded in an expanding universe, both for 
spatially-fiat and curved cosmologies, and compare our re- 
sults with McVittie's metrics and with work conducted by 
other authors on similar models. We also discuss the physical 
consequences of our derived metrics, focussing in particular 
on particle dynamics and the force required to keep a test 
particle at rest relative to the point mass. 

The outline of this paper is as follows. In Section 2, 
we introduce the tetrad-based method for solving the Ein- 
stein equations for spherically-symmetric systems, and de- 
rive the metrics for a point mass embedded in an expanding 
universe for spatially-fiat, open and closed cosmologies. In 
Section 3, we compare our metrics with those derived by 
McVittie. The geodesic equations for our spatially-fiat met- 
ric are derived and interpreted in Section 4. In Section 5, 



we derive a general invariant expression, valid for arbitrary 
spherically-symmetric spacetimes, for the force required to 
keep a test particle at rest relative to the central point mass, 
and consider the form of this force for our derived metrics. 
Our conclusions are presented in Section 6. 

We note that this paper is the first in a set of two. In our 
second paper (Nandra, Lasenby & Hobson 2011; hereinafter 
NLII2), we focus on some of the astrophysical consequences 
of this work. In particular, we investigate and interpret the 
zeros in our derived force expression for the constitutents of 
galaxies and galaxy clusters. 



2 METRIC FOR A POINT MASS IN AN 
EXPANDING UNIVERSE 

We derive the metric for a point mass in an expanding uni- 
verse using a tetrad-based approach in general relativity (see 
e.g. |Carroll|2003l ); our method is essentially a translation of 
that originally presented by [Lasenby et al.| ( |1998[ ) in the 
language of geometric algebra. First consider a Riemannian 
spacetime in which events are labelled with a set of coor- 
dinates x'^ , such that each point in spacetime has corre- 
sponding coordinate basis vectors e^, related to the met- 
ric via e.f_i ■ e.v = g^v At each point we may also define a 
local Lorentz frame by another set of orthogonal basis vec- 
tors e-i (Roman indices). These are not derived from any 
coordinate system and are related to the Minkowski metric 
riij — diag(l, —1, —1, —1) via • = rjij. One can describe 
a vector v at any point in terms of its components in either 



basis: for example «p = v-e^ and Vi = v-e.i. The relationship 
between the two sets of basis vectors is defined in terms of 
tetrads, or vierbeins eu^ , where the inverse is denoted e^ 
^ 

fip = e^^fifc. (6) 

It is not difficult to show that the metric elements are given 
in terms of the tetrads by gr^^ = Vij^^ ti^''v 

We now consider a spherically-symmetric system, in 
which case the tetrads may be defined in terms of four un- 
known functions /i(r, t), f2{r,t), gi{r,t) and g2{r,t). Note 
that dependencies on both r and t will often be suppressed 
in the equations presented below, whereas we will usually 
make explicit dependency on either r and t alone. We may 
take the non-zero tetrad components and their inverses to 



be 








60" = 






= - /252), 


61" = 


/2, 


e\ 


= -/2/(/l5l - f292). 


eo^ = 


92, 


1 

e 


= -92/iflgl - 1292), 


1 

ei = 


91, 


1 

e 1 


= - /25'2), 


2 

62 = 


1/r, 


2 

e 2 


= r. 




l/(rsin6l), 




— rsinO. 


In so 


doing, we 


have made 


use of the invariance 



(7) 



eral relativity under local rotations of the Lorentz frames to 
align £2 and £3 with the coordinate basis vectors £2 and £3 



at each point. It has been shown in Lasenby et al. ( 1998 1 
that a natural gauge choice is one in which f2 ~ 0, which 
we will assume from now on. This is called the 'Newtonian 
gauge' because it allows simple Newtonian interpretations. 
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as we shall see. Using the tetrads to calculate the metric 
coefficients leads to the line element 



9i 



92 



dt'^ + dr dt - 



-\dr^ -r'^dO.'^. (8) 
Si 



We now define the two linear difTerential operators 

Lt = fidt+ g2dr, 
Lr = gidr, 



(9) 



and additionally define the functions F{r,t), G{r,t) and 
M(r,t) by 



Ltgi = Gg2, 
Lrg2 = Fgi, 

M = \r (52^ - .gi^ + 1 - |Ar^ 



(10) 



where A is a constant. Assuming the matter is a perfect 
fluid with density p(r, t) and pressure p(r, t), Einstein's field 
equations and the Bianchi identities can be used to yield re- 
lationships between the unknown quantities, as listed below 
( Lasenby et al.|1998| ): 



Lrfl 

Lrgi 



Gfi fi 
M 



exp 



G 
5i 



dr 



Fg2 + 



|Ar — 47rrp, 



Ltg2 = Ggi 



M 



+ kAr — Anrp, 



Aiig2r p, 
232 



LtM - 

Ltp ■- 

LrM = 4:ngir''p, 
LrP = -G{p + p). 



+ F\{p + p), 



(11) 



From the LrM equation we now see that A'l plays the role 
of an intrinsic mass (or energy) interior to r, and from the 
LrP equation it also becomes clear that G is interpreted 
as a radial acceleration. In such a physical set up A is the 
cosmological constant. 

In order to determine specific forms for the above func- 
tions it is sensible to start with a definition of the mass 
M{r,t). For a static matter distribution the density is a 
function of r alone, p = p{r), and M{r) = Anr'^ p (r') dr' . 
Setting M(r) equal to a constant m leads specifically to 
the exterior Schwarzschild metric in 'physical' coordinates 
([3|. For a homogeneous background cosmology, p = p(t) 
and M{r,t) = (4/3)7rr^p(t), leading to the FRW metric in 
'physical' coordinates ([T]). In this work we choose M to de- 
scribe a point object with constant mass m, embedded in a 
background fluid with uniform but time-dependent spatial 
density: 



M{r,t) = ^■Kr-p{t) +m, 



(12) 



which is easily shown to be consistent with the LrM equa- 
tion above. 

We point out that the background fluid is in fact a to- 
tal 'effective' fluid, made up of two components: baryonic 
matter with ordinary gas pressure, and dark matter with an 
effective pressure that arises from the motions of dark mat- 
ter particles having undergone phase-mixing and relaxation 



(see Lynden-Bell (19671 and Binney fc Tremaine (2008 1 ) 



The degree of pressure support that the dark matter pro- 
vides depends on the degree of phase-mixing and relaxation 
that the dark matter particles have undergone, which (in the 
non-static, non-virialised case) will be a variable function of 
space and time. The properties of this single 'phenomeno- 
logical' fluid, with an overall density p{t), are studied in 
more detail in our companion paper NLH2. Here we simply 
calculate the total pressure p{r, t) of the background fluid 
required, in the presence of a point mass m, to solve the Ein- 
stein fleld equations in the spherically-symmetric case. The 
'boundary condition' on this pressure (at least for a flat or 
open universe, where spatial infinity can be reached) is that 
the pressure tends at infinity to the value appropriate for the 
type of cosmological fluid assumed. This is p(oo,t) = in 
the present case, since we are matching to a dust cosmology. 
The total pressure at finite r may be thought of as a sum 
of the baryonic gas pressure and dark matter pressure, but 
without an explicit non-linear multi-fluid treatment we do 
not break the fluid up into its components in the strong-fleld 
analysis presented below. 

Note that the central point mass in our model is in- 
evitably surrounded by an event horizon. The fluid con- 
tained within this region remains trapped and cannot take 
part in the universal background expansion, and so our ex- 
pression for M in ( |12[ ) is only valid outside the Schwarzschild 
radius. We thus expect the metric describing the spacetime 
to break down at this point. However, since one is usually 
most interested in the region m <^ r <^ R{t) (roughly equiv- 
alent to m ^ r ^ 1/H{t)), it is appropriate to continue 
using this deflnition for M to study particle dynamics far 
away from the central point mass. 



We are able to use ( 12 1 to determine specific forms for 



the tetrad components. We first substitute it into the LtM 



equation from (111 and simplify to obtain 



dp{t) 
^ dt 



3p2 



ip{t)+p). 



(13) 



Combining this result with the Ltp equation from (111 and 



the definition of F from ( 10 1, one quickly finds that 



g2 ^ dg2 
r dr 

This is easily solved for 32, and hence F, to yield 



52 - 
F ■ 



rH{t), 
H{t), 



(14) 



(15) 



where H{t) is some arbitrary function of t. Substituting 



these expressions into the Lrgi equation from (111, and us- 



ing the definition of M from ( 10 1 to fix the integration con- 
stant, one finds that 



2 1 
5i = 1 



2m 2 / N 



where we have defined the new function 



(16) 



(17) 



It should be noted that, by interpreting H{t) as the Hubble 
parameter, the three terms on the right-hand-side of (17 1 
correspond to —k/R^{t) via the Friedmann equation for a 
homogeneous and isotropic universe, where k is the curva- 
ture parameter and R(t) is the scale factor. Calculating the 
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function G is now straightforward from its definition in ( 10 1 



G = 



2H{t){r^r,{t)+m) 



(18) 



Finally, the function /i can then be calculated from the 



Lr/i equation in (11 1. We thus have expressions for all the 



required functions /i, gi, g2, F and G. 

We conclude our general discussion by noting the rela- 
tionship between the fluid pressure p and the function /i. 



Combining the Lr/i and LrP equations in (111, one quickly 
finds 



drfi drfi 

OrP = —^PW- 

Ji n 



(19) 



This first-order linear differential equation can be easily 
solved for p by finding the appropriate integrating factor, 
and one obtains 



(20) 



where ^ {t) is, in general, an arbitrary function of t. Combin- 
ing this result with (131, (15 I and (171, and recalling that 



r]{t) = —k/R (t), one quickly finds that 



1 

4tt 



dH{t) 
dt 



(21) 



Using the Friedmann acceleration equation for a homoge- 
neous and isotropic universe, one then finds that ^(t) — 
(1 -|- w)p{t), where w is the equation-of-state parameter of 



the cosmological fluid. Hence the relationship (201 between 
the fluid pressure and /i becomes simply 



p = p{t)[{l + w),h~l]. 



(22) 



For an FRW universe without a point mass (m — 0), in 
which /i = 1, we recover the relationship w = p/ p{t). 



2.1 Spatially- fiat universe 



A number of observational studies, such as WMAP (Larson 



et al. 20111, indicate that the universe is spatially flat, or 



at least very close to being so. In this case rj{t) — and 
the resulting expressions for the quantities g\, g2, F , G and 
/i are easily obtained; these are listed in the left-hand col- 
umn of Table [T] We note that the given expression for /i 
is obtained by imposing the boundary condition /i — >■ 1 as 
r — ^ oo; from \22\ this follows from the physically reason- 
able boundary condition that the fluid pressure p — >■ as 
r — >■ 00. 

From (|8| this leads to the metric (in 'physical', i.e. non- 
comoving coordinates): 



r 



dr + 2rH(t) 1 



2m. 



- 1- 



2m 



dr dt 
(23) 



which is a natural combination of ([T]) and ([3|. Indeed, it can 
be seen to tend correctly to the spatially-flat FRW solution 
([l| in the limit m — > (or r — >■ oo), and to the Schwarzschild 
solution (|3| in the hmit H{t) 0. 



In this case, the general expression ( 22 1 for the fluid 
pressure becomes 



Pit) 



2m 



(24) 



This can be checked directly by substituting our form for /i 
in the LrP equation from (111, from which it follows that p 
and p{t) are related by 



P + Pit) 



2m /r 



dr. 



(25) 



Imposing the boundary condition that the pressure tends to 
zero as r — > oo, this leads to (241, as expected. 

We note that the metric ( 23 1 is singular at r = 2m. 
This is, however, unlike the r = 2m coordinate singu- 
larity of the standard Schwarzschild metric. The latter 
arises due to a poor choice of coordinates and by convert- 
ing to another more suitable coordinate system, such as 
Eddington-Finkelstein coordinates, it can be shown that the 
Schwarzschild metric is actually globally valid. On the con- 



trary, for our derived metric, we see from ( 24 1 that the fluid 
pressure becomes infinite at r = 2m, which is thus a real 
physical singularity. Hence our metric is only valid in the 
region r > 2m. In reality, this region is usually deeply em- 
bedded within the object. In an attempt to make our so- 
lution globally valid, we shall present an extension of this 
'exterior' work to the interior of the object in a subsequent 
paper, where we will take its spatial extent into proper con- 
sideration and follow a similar approach to that used in this 



work. We point out that Nolan ( 1999b I has suggested con- 



sidering a different type of fluid altogether, such as a tachyon 
fluid, to deflne an equivalent metric inside this region, but 
we leave this type of approach for future research. 



2.2 Open universe 

For an open universe (k = —1), one has ri{t) = 1/B?{t) 
and the resulting expressions for the quantities g\, g2 and 
F are easily obtained and are listed in the middle column 
of Table [l] In this case, however, the expression for /i (and 
hence G) is less straightfoward to obtain. Combining ( 18 1 
with the Lrfi equation from ( |11[ |, one finds that /i may be 
written analytically in terms of an elliptic integral: 



1 



R^it) 



2m 



RHt) 



r dr 



W{t) 



3/2 ■' 



(26) 



where the constant of integration, or equivalently the limits 
of integration, must be found by imposing an appropriate 
boundary condition. 

To avoid the complexity of elliptic functions, we instead 
expand /i as a power series in m, since astrophysically one 
is most interested in the region m <^r <^ Rit), i.e. values of 
r lying between (but far away from) the central point mass 
and the curvature scale of the universe. Recasting /i in its 
differential form gives 
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Table 1. Functions defining the metric for a point mass embedded in an expanding universe for a flat (fe = 0), open {k = —1) and closed 
(fc = 1) cosmology. 



G 



k = 



k = -1 



k = l 



92 rHit) 
F H(t) 

% 1 



r ij2(t) 

H(,t) 



/-I 2m 7-^ 
V r R2(t) 

rH(t) 

H(t) 




200 400 600 



1000 1200 1400 1600 

r 



Figure 1. The logarithm of the fractional error in the approxi- 
mate series solution ( |28| for /i, relative to the exact numerical 
result ( |27| , for fixed arbitrary values in = 1 and R{t) = 1000. 
Numerical precision effects are visible beyond r Ri 1500. 



1 dfi ^ ^ m _ 2m ^ r 



/i 9r 



(27) 



The series solution to this equation is 

/i = 1 + ^ + 1^ + mniVRHt)+r^ + 0{m'), 

where the arbitrary function /3(t) can only be determined 
by the imposition of a boundary condition. For an open 
universe we expect p — >■ and hence G — >■ as r — > oo. 



Therefore, from (221 we expect /i — >■ 1 as r — > oo. This 
gives I3{t) = and hence 



, , , , m 2mr 
Mr,t) = 1+- + 



2m 



r R^{t) R{t) 



''' +0(m^). (28) 



R^t) 



As shown in Fig.[T] expanding to first order in m is sufficient 
to represent the solution to high accuracy for our region of 
interest. 

An approximate form for the metric in the case of an 
open cosmology is therefore 

ds^ = (joo dt^ + 2goi dt dr + gn dr^ — r^dil^ , 



where 
Poo « 

501 ~ 

511 = 



+ 



s^w) ( 



1 + ? + 



B2(t) 



2m 
R(t) 



1 + 



-1 2m I 



rH{t) 



2m r 



1 + - + 

r 

2 \ -1 



r_ _ 2m / 
t) R(t) Y 



1 + 



ii2(t) 

(29) 



It can be verified that in the limit m — >■ (or r — >■ oo) 
this reduces to the standard k = —1 FRW metric. Also, in 
the limit r/R{t) — >■ and working to first-order in m/r, the 



case ( 23 \ 



metric coefficients in 1 29 1 reduce to those in the spatially-fiat 



We also note that the metric is not singular at r = 2m, 
but instead becomes singular where 



2m 



+ 



- 0. 



(30) 



Indeed, fi is singular there. Multiplying through by r, the re- 
sulting cubic equation has a positive discriminant and hence 
only one real root, which occurs at a radial coordinate in- 
side the standard Schwarzschild radius r — 2m. Since /i and 
hence the fiuid pressure are singular there, then, as in the 
spatially-fiat case, this is a true physical singularity rather 
than merely a coordinate singularity. We further point out 
that, in contrast to the spatially-flat case, the radial coordi- 
nate at which this singularity occurs is a function of cosmic 
time t. 



2.3 Closed universe 

For a closed universe {k — 1), one has ri{t) = —1/R?{t) 
and the resulting expressions for gi, 32 and F are listed in 
the right-hand column of Table [T] As in the open case, the 
expression for /i (and hence G) requires more work. One 
finds that /i can similarly be given analytically in terms of 
an elliptic integral: 



1 



R^[t) 



2m 
r 



r dr 



1 - 



3/2 ■ 



(31) 



where, once again, the constant of integration or limits of in- 
tegration, must be found by imposing an appropriate bound- 
ary condition. As we will see below, however, the imposition 
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of such a boundary condition requires considerable care in 
this case, since the hmit r — >■ oo is not defined for a closed 
cosmology. Recasting /i in its differential form gives 

It is again convenient to expand /i as a power series in m, 
which reads 

/i = 1 + ^ - 1^ + mWR^t) - r2 + O(m^), (33) 

where the arbitrary function /3(t) can only be determined 
by the imposition of a boundary condition, to which we now 
turn. 



2.3.1 Boundary condition on /i 

The main problem in defining an appropriate boundary con- 
dition for a closed universe is that our 'physical' radial co- 
ordinate r only covers part of each spatial hypersurface at 
constant cosmic time t. One can see from Table [T] that gi 
and hence the metric becomes singular when 



2m 



= 0. 



(34) 



This also corresponds to where /i becomes singular, from 



equation (311, and hence where the pressure diverges, from 
equation (221. Multiplying through by r, the resulting cubic 
equation has a negative discriminant and hence three real 
roots, provided m ^ R{t) / (?>^/?i) . It is easily shown that one 
of these roots lies at negative r, and is hence unphysical, 
and the remaining two roots lie at 



n(t) 
r2{t) 



2R{t) 

2R{t) 
^/3 



1 -i/3\/3 m 

- cos , 

3 V Ri't) 

1 _i / 3\/3m 

- cos , 

3 V Ri't) 







+ 










+ 


6 



(35) 



It is straightforward to show that r\{i) corresponds to the 
'black-hole' radius, and lies outside the Schwarzschild radius 
r = 2m. At this point /i, and thus the fluid pressure, are 
also singular, and so this corresponds to a true physical sin- 
gularity, as in the spatially-flat and open cases. 

The other root, r2{t), is easily shown to correspond to 
the 'cosmological' radius, and lies insidethe curvature radius 
r — R{t). This radius should be merely a coordinate singu- 
larity, which we verify in Section [2. 3. 2| Thus one would not 
expect the fluid pressure, and hence /i, to be singular there. 
Moreover, one would also expect dfi/dr to be non-singular 
there. Prom the expression ( |32[ ), one quickly finds that for 
the latter condition to hold, one requires 



Mr2,t) 



r2{t) — 3m 
r2{t) - 2m' 



(36) 



Thus, the series expansion of /i about the cosmological ra- 
dius takes the form 



h{r,t) 



r2 



3m 



r2 



2m 



^a„(f)(r2-r)^ 

n = l 



(37) 



where the coefficients a„(t) may be determined by substitu- 
tion into ( 32 1, and we have momentarily dropped the explicit 



dependence of r2{t) and R(t) on t for brevity. One finds that 
the first coefficient, which is the only one of interest, reads 



ai{t) 



2 I r2-3m 

^rnr2 tt- 

r2 — 2m 



2(m_R2 -|_ - r2R^){r2 - 2m) - r^{r2 - 3m) ' 



(38) 

Thus, we have determined the (finite) values of both /i and 
dfi /dr at the cosmological radius r = r2. Since the differen- 
tial equation ( |32[ | for /i is first-order in its radial derivative, 
one can thus, in principle (or numerically), 'propagate' /i 
out of the cosmological radius, towards smaller r values. 
The boundary conditions ( |36[ ) and ( |38[ ) therefore uniquely 
determine /i. 

We point out that, in addition to being singular at the 
'black-hole' radius r = ri{t), the function /i (and hence 
the fluid pressure) will also be singular at the zeros of the 
integral given in equation (31 1. If the inner-most zero occurs 
at r, (t), which is some (unique) function only of m and R{t), 
we may thus represent /i in the integral form 



1 



2m 



u du 



r,(t) M _ 2m _ 



3/2 ■ 



(39) 

It is not clear how to find an analytical expression for r,{t), 
but numerical results show that r = r, (t) lies slightly outside 
the radius r = ri(t). Moreover, as m/R(t) — > 0, both the 
absolute and fractional radial coordinate distance between 
the two radii decreases. Indeed, in any practical case, the 
two will be indistinguishable. 



Turning to the power series approximation ( 33 1 of /i , 
valid for our region of interest m <^ r <^ Rif), the identi- 
fication of the appropriate boundary conditions at the cos- 
mological radius r — r2{t) now allows us to fix the arbitrary 
function P{t) straightforwardly. From the small m approxi- 
mation of ( 34 1 , it is clear that the limit r — > r2 (t) is equiva- 



lent to r approaching R{t) from below. If /3(i) were non-zero, 
then dfi/dr would be singular at the cosmological radius, 
owing to the ^/ R'^{t) — 
require /3(t) =0, so that 



term. We thus deduce that we 



/i 



1+^ 
r 



'irm 2 s 

+ 0(m 



RHt) 



(40) 



As shown in Fig. [2] this expansion is sufficient to represent 
the solution to high accuracy for our region of interest. Thus, 
an approximate form for the metric in the case of a closed 
cosmology is 



ds^ = goo dt^ + 2goi dt dr + gn dr^ 



where 



500 



r^H^t) 



1 - 



501 



5ii 



rH{t) 



1 I m 2mr 



W(t) 



^ + T ^ Rmt) 



2m 



(41) 



r R^(t)^ 

It can be verified that in the limit m — > 0, this reduces to the 
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Figure 2. The logarithm of the absolute fractional error in 
the approximate series solution | |40[ l for /i, relative to the ex- 
act numerical result \39\ , for fixed arbitrary values m = 1 and 
R{t) = 1000. The spike at r 580 is the result of the fractional 
error changing sign there. 



standard fc = 1 FRW metric. Also, in the limit r / R{t) — s> 
and working to first-order in m/r, the metric coefficients in 
(411 reduce to those in the spatially-flat case (23 1. 



2.3.2 Alternative radial coordinate 

We have mentioned that, for a closed universe, our 'physical' 
radial coordinate r only covers part of each spatial hypersur- 
face at constant cosmic time t. The metric has a singularity 
at the 'cosmological' radius r — r2{t), which we now verify 
is merely a coordinate singularity. 

Even without the presence of a point mass, a similar 
problem arises when using the 'physical' r-coordinate in the 
case of a pure closed FRW metric, which has a coordinate 
singularity at r = R{t). This issue is discussed in |Lasenby| 



et al. ( 1998), where an alternative radial coordinate was in- 



troduced, which removes the singularity at the cosmologi- 
cal radius. The solution presented there amounts to a two- 
stage coordinate transformation: one first transforms to a co- 
moving radial coordinate and then performs a stereographic 
transformation. The final form of the metric is 'isotropic' in 
the sense that its spatial part is in conformal form. Thus, 
an obvious approach in our case (with m 7^ 0) is to seek 
an isotropic form for the metric, which reduces to the form 
found by Lasenby et al. ( 1998 1 when m — 0. 

We note that only the radial coordinate r is trans- 
formed, so the t coordinate keeps its meaning as cosmic time. 
We thus consider a new radial coordinate of the general form 
f = r{t,r). In fact, it is more convenient in what follows to 
consider the inverse transformation r = r(t,f). It should be 
understood here that r is a new function to be determined, 
the value of which is equal to the old radial coordinate. 

We begin by considering the metric in the form (|8|, 
where the functions g\ and g2 are given by the analytical 
expressions given in the right-hand column of Table [l] For 
the moment, we will not assume a form for /i. Performing 
the transformation r = r{t,r), we will obtain a new metric 
in t and f (and the standard angular coordinates 9 and 0). 



By analogy with the approach of Lasenby et al. (19981, we 
require this new metric to be in isotropic form, i.e. the co- 
efficient of df^ must equal that of dd^ and the cross-term 



dtdr must disappear. The first condition leads to 
dry _ r[rR^{t) -r^ -2mR^{t)] 

df-j ^ -prnt) ' 

while the second condition yields the following direct for- 
mula for /i: 

fi=' 



1 dlnr 
H{t) dt 



(42) 



(43) 



When these conditions are satisfied, the resulting metric is 
given by 



ds^ 



f 1 91nr 

\mt) dt 



dt^ 



— [df^ +P(dr -hsin 



(44) 

which depends only on the single function r = r{t,r). Sub- 
stituting this metric into the Einstein field equations, one 
can verify that it is consistent with a satisfactory perfect 
fiuid energy-momentum tensor, in which the fluid density 
depends only on t and the radial and transverse pressures 



are equal and satisfy the relation ( 22 1 



One can solve ( 42 1 to obtain an expression in integral 
form for r in terms of r. There are two solutions, one of 
which reads 



Inr 



R{t) du 



(45) 



r-At) \/u\uR^{t) -u^ ~ 2mR{t)] 

Here we have used r, (t) as the lower limit of integration, so 
that we can more easily make a connection with our inte- 



gral expression for /i given in ( 39 1 . This might not be the 



appropriate integration limit in this case, however, and so 
we include the arbitrary function C,{t) to absorb any discrep- 
ancy. 

The arrangement of integration limits in ( |45| l enables 
us to consider the r-range from near the point mass to the 
cosmological radius. One can now proceed beyond the cos- 
mological radius, however, by writing the second solution to 
(421) as 



Inr 



■At) 



R{t) du 



(t) Jr,(t)) ^Ju[uR^{t) - m3 - 2mR{t)] 



+ CW, 
(46) 



which also satisfies (42) and reduces to (45) at r = r2{t). In 



Fig. [3] we plot Inr over the full range of r. The maximum 
value of r occurs at the cosmological radius and, thereafter, 
f continues to increase as r decreases again. The geometrical 
interpretation of this result is illustrated in Fig. [4] with one 
spatial dimension suppressed. In essence, as f increases, one 
is following a great circle on the surface of a sphere, starting 
at the original point mass and ending at an image mass at 
the antipodal point of the universe. 

We note that there is an interesting relationship be- 
tween any two r values that correspond to the same value 
of r. If two such values are fa and ft, then 



ln{fafb) 



2(t) 



R{t) du 



-2C(t). (47) 



r,(t) y/u[uR^{t) -u3 - 2mR{t)] 

Since the right-hand side is a function only of t, then at 
any given cosmic time the values fa and ft are reciprocally 
related. This behaviour also occurs in the case of a pure 



closed FRW model, as discussed in Lasenby et al. ( 1998 ) 
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Figure 3. The relationship l |46| between the radial coordinates r 
and f for a closed universe in arbitrary units (see text for details). 



the latter), we can evaluate the above integral for C^(t). The 
only remaining ambiguity is an arbitrary additive constant, 
which corresponds to an arbitrary multiplicative constant 
in the definition of f in terms of r, i.e. it is not possible to 
identify a unique overall scale for r, which seems reasonable. 

We note that our solution involving an 'image mass' 
at the antipodal point of the universe ties in with the sce- 



nario recently investigated by Uzan, Ellis & Larena (20111, 



who also considered a closed universe with masses embed- 
ded symmetrically at opposite poles. They showed a static 
solution was not possible for this case, which fits in well 
with the fact that here we have an explicit exact solution 
for the masses embedded in an expanding universe. The ex- 
act nature of the correspondence with their work will be the 
subject of future investigation, however. 




Figure 4. Illustration (with one spatial dimension suppressed) of 
the geometry of a spatial hyperface in the closed universe solution. 
The solid concentric circles around the original point mass m 
represent, respectively, the 'black-hole' singularity at r = ri(t) 
and the surface r = rt,{t), at which the fluid pressure diverges. 
The circle around the equator represents the cosmological radius 
r = r2(t), in which the geometry is 'reflected' to yield an image 
mass m at the antipodal point in the universe. 



So far, we have left (^{t) undetermined, but one can in 
fact obtain an expression for ({t) by combining the standard 
partial derivative reciprocity relation 



dr L 



with the expression ( 43 1 for fi , which gives 



h = ~rH{t) 



dr 



(48) 



(49) 



Finding the derivatives of r from ( 45 1 and equating this 



result with our integral expression ( 39 1 for /i then yields 



at) 



R{t) 



dr. 



^r,[R^{t)r, -ri- 2mR'^t)] dt 



dt + constant. 



(50) 



Although we do not have an explicit expression for rt{t), 
we have shown that there is an operational method for de- 
termining it (i.e. where /i becomes singular when numeri- 
cally propagating it inwards from the cosmological radius). 
If we know r, (t) and R{t) at each time slice (see below for 



2.3.3 Cosmological evolution 

When considered as a function of the new radial coordinate 
f, the radial derivative of /i is given by 



dr 



dr 



(51) 



In Section 2.3.1 we showed that dfi/dr is finite at the 
cosmological radius r = r2{t), whereas (421 shows that 



dr/df — there. Thus, we conclude that dfi/dr = at this 
point. From (221, this corresponds physically to the gradi- 



ent dp/df of the fluid pressure vanishing at the cosmological 
radius. It appears, therefore, that the fluid pressure at this 
point can be any function of t alone; once this is specified, 
one has sufficient information to solve jointly for /i and R{t). 
This seems plausible physically, since we have not specified 
an equation of state relation between p and p, and so need 
instead to specify a boundary condition on p at the cosmo- 
logical radius. 

Let us consider the specific case in which we impose as 
our boundary condition p{r,t) = at r = f2(t), given by 
([35|, for all t. Remembering that fi{r2,t) is given by ([36| 
and using ( 22 1 and the standard cosmological field equa- 
tions, one can then show that 



R"{t) = 



Rit) 



1 + R'\t) 



4 cos^ cos~^ 



R(t) 



(52) 



We can therefore, in principle, obtain the expansion history 
R{t) by solving this second-order differential equation. In 
the limit m/R{t) <C 1, an approximate first-integral of the 
equation is given by 



R!^{t) 



C / 3m \ 



(53) 



where C is a constant. Substituting this result into the ex- 



pression for p given in ( 57 1, the corresponding expression for 



the fluid density is given by 



3C / 3m \ 

9,11 R?(t) """^^ yW)) 



(54) 



This shows that the presence of the point mass m means 
that p{t) does not quite dilute by the usual 1/R?'{t) factor. 
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Table 2. A comparison of the fluid pressure and density obtained 
from McVittie's metric and our metric for a point mass m in a 
spatially-flat expanding universe (with A = 0). 



McVittie 



This work 



87rp 
87rp 



-3H^(t) - 2H'{t) 



2H'{t) 



3 COMPARISON WITH MCVITTIE'S METRIC 

At the same time that McVittie derived his metric ([5| for a 
mass particle in a spatially-flat expanding universe, he also 
attempted to extend his result to apply to a universe with 



arbitrary spatial curvature ( McVittie|1933 19561 



-'- 2rfl(t) V 1 ' 4 



1+^ 



2rfl(t) 



dr 



2rR(t) 



1 + ^ 



-i?^(t)(dr^ (55) 



We now compare our metrics for k = 0, k = —l and k = 1 
with this result. 



3.1 Spatially-flat universe 



We first note that our A: = metric ( 23 1 behaves similarly to 
McVittie's metric ([5| in the appropriate limits, as it should. 
In fact, we now show that our metric and McVittie's are re- 
lated by a coordinate transformation, despite the fact that 
they were derived in very different ways. We deduce this re- 
lationship by comparing specific physical quantities. Assum- 
ing a perfect fluid energy- momentum tensor, precise forms 
for the density and pressure can be derived from the metrics 
(Carrera & Giulini 2010a b), without assuming the relation- 
ship 1 22 1, which may not hold for McVittie's metrics. The 



resulting quantities are given in Table [2] assuming for sim- 
plicity that A = (although our conclusions still hold in the 
A 7^ case). 

It is found that the background densities already match 
in both models and correspond to the density obtained from 
the FRW metric; this is not suprising since we are still work- 
ing within an FRW background. We can use the different 
forms for the pressure to define a coordinate transforma- 
tion: 



r = fR{t) 1 + 



2rR{t) J 



(56) 



It is easily shown that our metric is equal to McVittie's 
metric under this transformation. The transformation actu- 
ally converts our 'physical' radial coordinate to its comov- 
ing and isotropic analogue in a single step. Hence our met- 
ric describes the same spacetime as McVittie's metric, with 
the bonus that its interpretation is more transparent in our 
non-comoving coordinates. Note that the limit r — > oo corre- 
sponds to r — >■ oo, leading back to the FRW pressure in both 



cases, as expected. The limit r — >■ on the other hand, which 
also corresponds to r — >■ oo, is not well-defined. As pointed 
out by|Nolan|(|1998||1999a|b| ), and more recently [Faraoni fc 
Jacques ( 2007| ), it is clear that there is a spacelike singular- 
ity in McVittie's metric at r = m / (2R{t)) , corresponding to 
a diverging pressure. The interpretation of this singularity 
has been under debate for some time, but we now see that 
it corresponds simply to the physical singularity at r = 2m 
in our coordinate system. This in turn coincides with the 
location of the event horizon from which the background 
fluid is unable to escape, as discussed in Section |2.1[ Thus 
McVittie's metric can only be valid for r > m/{2R{t)). 

This transformation to 'physical' coordinates has also 
previously been pointed ou t by |Nolan| |T998| |1999a| and 



other authors ( Arakida 2009 Bolen, Bombelli & Puzio 



2001 



Faraoni & Jacques 2007). They have also highlighted the 'ac- 



cident' by which McVittie initially derived his metric, which 
in its original form does not describe a central mass. The 
mass m is located at r = in our setup, but this does not 
correspond to the radial coordinate f — 0. Instead the point 
mass is located at r = —m/{2R{t)), which seems rather un- 
natural. McVittie's accident led to problems in the cases of 
spatially curved cosmologies, as we shall see shortly. 

used the transformation ( 56 1 to allow 



Nolan (1999a 



for a more intuitive analysis of the global properties of 
McVittie's metric. In particular, Nolan identified the func- 



tion M{r,t) in ( 12 1 as the Misner-Sharp energy of the space- 
time ( Misner fc Sharp|1964 l. This is a measure of the 'total 
energy' of each fluid sphere in terms of the work done on it 
by the surrounding fluid. We note that this was the start- 



ing point in our derivation of the metric (231, as opposed to 



an emergent feature. It therefore clarifies a point that has 
been under debate for a long time; McVittie's metric does 
indeed describe a point-mass in an otherwise spatially-fiat 



and Lake fc Abdelqader (20111) 



FRW universe (see also Kaloper, Kleban fc Martin (20101 



3.2 Spatially-curved universe 

We have already shown that an advantage of our approach 
is that we have a natural way of generalising to the case of a 
curved cosmology. For = ±1, our model still incorporates 
an FRW background so we would again expect our form 
for p{t) not to deviate from the standard FRW result; in 
particular, it should be a function of t only. Indeed, using 
the general metric ([8| with gi and g2 given in Table [l] we 
find (again assuming A = for simplicity) 



8t: p{t) = 3 {H\t)+ ^ 



8np = -3 H\t) + 



R^{t) 
k 



RHt) 



-2[H'{t) 



RHt) 



where the series solutions for /i for the cases k 
k = 1 are given in ( 28 1 and ( 40 1 , respectively. 



(57) 



-1 and 



For McVittie's metric (551, however, we find the pecu- 



liar feature that the background density does depend on the 
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radial coordinate r: 



87rp = 3ff^(t) + 



3A; 

mt) 



1 + 



^■Kp = -ZH'^it) -2H(t) 



m \ 
2wR{t) ) 

1 H T: ■ 

' 2wR(t) 

-I m 

2wR(t) 



where w = 



(1 + 



1/2 



(58) 



(Nolan 19981. This is odd since 



t) ) (l + 2™H(t) 



the only apparent difference between this model and the 
spatially-flat case is the spacetime curvature. McVittie's 
fc = ±1 metrics hence cannot describe a mass particle em- 
bedded in a background with homogeneous density. Indeed 
the density and pressure ( 58 1 do not even asymptotically 
tend to the FRW solutions. 

Our fc = ±1 metrics, which are derived by explicitly as- 
suming a model of a point mass embedded in a homogeneous 
background, are therefore inherently different to McVittie's 
solutions. |Nolan| (jl998 ) has used a model similar to ours to 
obtain a metric in the fc = — 1 case that is expressed in terms 
of an elliptic function, but we have not yet found a coordi- 
nate transformation that equates this with our solution. A 
full comparison will be the subject of future research. 



4 GEODESIC MOTION IN THE 
SPATIALLY-FLAT METRIC 

The motion of a test particle moving under gravity around 
a central point mass in a static spacetime has been very well 
studied. In simple terms, we would expect the incorporation 
of an expanding background to provide an additional 'force' 
that alters the trajectory of the test particle. We investigate 
this possibility by first calculating the geodesic equations 
for our metric, using the usual 'Lagrangian' technique. We 
restrict our attention to the fc = case, since observations 



indicate the universe to be very close to spatially flat ( Larson 
et al.|2011| |. 

Working in the equatorial plane 6 = n/2, the 'La- 
grangian' corresponding to our flat metric ( 23 1 is 



1 - 



2m 



2rH{t) 

\ _ 2m 



where dots denote differentiation with respect to the proper 
time r of the test particle. The three remaining geodesic 
equations are obtained from the Euler-Lagrange equations 
— {§iii) for x'^ — t, r and 0. For stationary, 
spherically-symmetric metrics, such as the Schwarzschild 
metric, one can immediately obtain the first integrals of the 
t and (j) geodesic equations. Moreover, it is usual to replace 
the r geodesic equation by its first integral Qii^x'^i" = 1 (for 
timelike geodesies) . These three equations can then be easily 
combined to obtain an 'energy equation' in terms of r, r and 
constants of the motion only. By differentiating this equa- 
tion, a simple expression for r can be obtained, if desired. 
This procedure is followed for the Schwarzschild de Sit- 
ter metric by |Balaguera-Antoh'nez, Bohmer fc Nowakowski| 
(20061. Our metric, however, is a function of both r and t, 



and therefore this standard approach does not prove useful 
in finding r and the procedure must be modifled. 



In our case, the r geodesic equation is not replaced since 
it contains an r term. This is of particular interest to us in 
determining the nature of the gravitational 'forces' acting 
on the particle. The geodesic equation is r^^ = L, where 
L is the speciflc angular momentum of the test particle. 
Substituting this into the r equation and rearranging gives 



2m 
r J r 

2m 

r 



2m 



rH'{t)P + 



rH^{t)]r + rH{t)Jl 



2m. 



(59) 



The t term can be eliminated completely using the t- 



equation, so that (59 1 now becomes 



2m 



r'H\t) 



2m 
r 



+ 



1 



- 2rH{t)- 



rt. 



(60) 



Note that the i terms above, when combined, differ very 
slightly from that derived in |Carrera fc Giulini| ( [2010b[ ); we 
believe the latter to be in error. We now use the condition 
gfiuX^x'^ = 1 to eliminate the terms. After some algebraic 
manipulation this finally leads to a relatively simple exact 
expression for r, namely 

( ^ 3m \ m ,,2/ s 

^=7T -I2+-^'W 



2r^H{i)H'{t) 



1 - 



rt 



+ 



+ 



1 - 



2m 



rH'it) 



1 + 



2m 



2m 



(61) 



4.1 Newtonian limit and forces 



It is of interest to consider special cases of (611. First we 



consider the weak field approximation by assuming m ^ 
r <^ 1/H{t) and expanding to leading order in small quan- 
tities. We also take the low-velocity limit t ~ 1, r ~ 0. In 
this Newtonian limit, we find 



dt^ ^ ^ ' ^ ' J.2 J.3 

R"(t) m 



R(t) 



(62) 



where q(t) = —H'{t)/H^{t) — 1 is the deceleration pa- 
rameter. This result can also be obtained directly from 
our fiat metric (|23[) in its Newtonian limit ( |Nesseris fc| 
Perivolaropoulos|2004 1 : 



ds 



2m 
r 



dt^ + 2rH{t)drdt - dr^ - r'^d^f. 

(63) 



Note that this incorporates the usual low velocity approxi- 
mation, but the weak field condition used is simply m/r <C 
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1. Incorporating the condition rH{i) ^ 1 as well simply 
leads to the flat Minkowski metric, which does not effectively 
describe the physical system in which we are interested. 



We may interpret ( 62 1 as the physical acceleration of 



the test particle in the Newtonian limit, and hence the grav- 
itational force per unit mass acting upon it. Aside from the 
'centrifugal' term depending on the specific angular momen- 



tum L, the terms in (621 correspond to the standard l/r 



inwards force due to the central mass and a cosmological 
force proportional to r that is directed outwards (inwards) 
when the expansion of the universe is accelerating (deceler- 
ating) . 

An interesting feature of the cosmological force is that 
its direction (and magnitude) depends on whether the uni- 
verse is accelerating or decelerating (as determined by the 
sign of q{t)). This was pointed out previously by Davis, 



Lineweaver & Webb ( 2003 1 , and highlights the common mis- 



conception of there being some force or drag associated sim- 
ply with the expansion of space; instead the force is better 
associated with the acceleration/ deceleration of the expan- 
sion. In particular, we note that if the universe is decelerat- 
ing, the cosmological force is directed inwards and the test 
particle inevitably moves towards the central mass m, falls 
through its position and joins onto the Hubble flow on the 
other side. This does not, however, argue against the idea of 
the expansion of space. 

Since the current deceleration parameter go of the uni- 
verse is measured to be approximately —0.55 (Larson et al 



2011), we see that the present-day cosmological force is di- 



rected outwards, as one might have naively expected intu- 
itively. However, since the expansion history of the current 
concordance model of cosmology changes from a deceler- 
ating phase to an accelerating phase at a relatively recent 
cosmic time, the cosmological force must reverse direction 
at this epoch. A realistic scale factor corresponding to the 
standard spatially-flat concordance model is ( [Hobson, Efs-| 
tathiou fc Lasenby|2006 1: 



(l-»A,o) 



sinh^ (^^Ho\/^A,ot 



m 

Ro 

where Qa^ ~ 0.7 ( Larson et al.|201l| ) is the current fraction 
of the critical density of the universe in the form of dark 
energy. Using this expression, in Fig. [5] we plot the ratio of 
the cosmological force relative to its present-day value, as a 
function of redshift z. As anticipated, we see that the cos- 
mological force reverses direction; it is an inwards force for 
redshifts larger than about z = 0.67. Moreover, the magni- 
tude of the force relative to its current value increases by a 
factor of 10 by about 2; = 2.5 and by a factor of nearly 60 by 
z — 5; it continues to grow with redshift quickly thereafter. 

The cosmological force appears to have been cor- 
rectly taken into account in cosmological simulation codes. 



For example, in Springel, Yoshida & White (20011 and 
[Springei] ( |2005[ ), dark matter and stars are modelled as 
self-gravitating coUisionless fluids following the coUisionless 
Boltzmann equation, and this term appears in the equation 
of motion of the particles. For a single central mass m, the 
equation of motion of a test particle is taken as 



m 



r2i?3(t)' 

where f represents comoving coordinates. 



(64) 

Converting to 




Figure 5. The ratio of the cosmological force at redshift z to its 
current value. 



'physical' coordinates via r — fR{t) recovers the result ( 62 1 
with i = 0. 

We perform a full exploration of the astrophysical con- 
sequences of the geodesic equation ([62|) in our companion 
paper NLH2. We also note thatlPricc (2005) has used a very 
similar geodesic equation to analyse the motion of an elec- 
tron orbiting a central nucleus, with a slight modification to 
account for the dominant electrostatic force in that problem. 
An extension of this work is also presented in NLH2. 



4.2 Particle released from rest 

It is, of course, of interest to calculate the 'force' experienced 
by the particle in the fully general-relativistic case. This is 
discussed in detail in the next section, but as a prelude let us 



return briefiy to (61) and consider a particle released from 
rest, for which r = L = Q. Thus, instantaneously, one has 



r = rH^{t) 



m 



rH\t){q{t) + l)^ 
1 - ^ - r2_ff2(t) 



(65) 



One should bear in mind at this point, however, that the 
'physical' coordinate r is not a proper radial distance. Hence 



one cannot simply interpret ( 65 1 as the instantaneous force 



on a particle released from rest, or equivalently the negative 
of the force required to keep a test particle at rest relative to 
the central point mass. The proper expression for this force 
is derived in the next section. 



5 FRAMES AND FORCES 

We now derive a fully general-relativistic invariant expres- 
sion, valid for any spherically-symmetric system, for the 
radial force required to hold a test particle at rest rel- 
ative to the central point mass. We first point out that 
the relationship between the 'physical' coordinate r and 
proper radial distance I to the point mass is defined through 
dl — —{l/gi)dr. It is therefore only in the spatially-fiat case, 
for which gi — gi{r) is independent of t, that r = or 
£ = are equivalent conditions. In general, one must thus 
choose which condition defines 'at rest'. Here we will adopt 
the condition r — 0, which corresponds physically to keep- 
ing the test particle on the surface of a sphere with proper 
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area iirr^. In practice, it would probably be easier for an 
astronaut (i.e. test particle) to make a local measurement 
to determine the proper area of the sphere on which he is 
located, rather than to determine his proper distance to the 
point mass, particularly in the presence of horizons. With 
this proviso, our derivation is valid for arbitrary tetrad com- 
ponents /i, gi and 32 (but we assume the Newtonian gauge 

/2 = 0). 

In general, the invariant force (per unit mass) acting on 
a particle is equal simply to its proper acceleration, i.e. the 
acceleration of the particle as measured in its instantaneous 
rest frame. We must therefore find an expression for the 
proper acceleration of a particle held at rest relative to the 
central point mass. This may be calculated directly in the 
coordinate basis, but it is more instructive for our purposes, 
and simpler, to perform the calculation in the tetrad frame. 

5.1 Tetrad frames and observers 

The tetrad frame defines a family of ideal observers, such 
that the integral curves of the timelike unit vector field eo 
are the worldlines of these observers, and at each event along 
a given worldline, the three spacelike unit vector fields Si 
[i = 1, 2, 3) specify the spatial triad carried by the observer. 
The triad may be thought of as defining the spatial coor- 
dinate axes of a local laboratory frame, which is valid very 
near the observer's worldline. In general, the worldlines of 
these observers need not be timelike geodesies, and hence 
the observers may be accelerating. 

Let us first consider the four-velocity and proper ac- 
celeration of the observer defined by our tetrad frame ([7|. 
The four- velocity of the observer is simply n = eo, so that, 
by construction, the components of the four-velocity in the 
tetrad frame are [ff] — [1,0,0,0]. Since — e^^u^, the 
four-velocity may be written in terms of the tetrad compo- 
nents and the coordinate basis vectors as u — fieg + g2ei. 
Thus, the components of the observer's four-velocity in the 
coordinate basis are simply [u^] — [t,r,6,(j)] — [/i, 32, 0, 0], 
where dots denote differentiation with respect to the ob- 
server's proper time. Hence, for all our previously derived 
metrics, we see from Table [l] that r = rH{t). Since the co- 
moving radial coordinate is defined through f = r/R(t), we 
deduce that f = rH{t){l — f\)/R{t). Hence the observer is 
not comoving with the Hubble fiow. Instead, since 1 — /i 
is negative, the observer's comoving radial coordinate is de- 
creasing. 

This behaviour is due to the presence of the cen- 
tral point mass, which results in our observer not moving 
geodesically (except as r — >■ 00). This is easily seen by cal- 
culating the proper acceleration a of our observer, which is 
given by a = a ■ a, where a = li is the four-acceleration 
of the observer. 

It is straightforward to show that, for a body moving 
with general four- velocity v, the four-acceleration is given 
in terms of the coordinate basis and the tetrad basis respec- 
tively by 



are known as Ricci's coefficients of rotation or the spin- 



(66) 



where V^^a are the connection coefficients corresponding to 
the metric ([sjl and 



connection. It has been shown by, for example, Kibble ( 1961 \ 



that these can be written in terms of the tetrad components 

as 



k fJ. V / r\ k f) ^ \ 



(68) 



It follows that ujijk is anti-symmetric in the first two indices, 
and Cijk is anti-symmetric in the last two indices. For general 
radial motion v'^ — = 0, in which case the components of 
the four-acceleration in the tetrad frame reduce to 

a — V +UJ iqV V + ii(f ) , 

-1 -1,1 /-0\2 I 1 -0-1 /en^ 

a — V +UJ oo(t^ ) + <^ oiV V , (69) 

and = = 0. It follows that for general radial mo- 
tion only cjioo = cooi (= -i^^o = ^i^^oo) and a;oii = Cno 



(= Lo II — u! Qi) need to be computed. Using (68 1, the defi- 
nitions of the tetrads and their inverses ([7| and the relations 
given in equation (llOl), one can show that 



c^oi ~ ^drfi ~ — G => a;ioo = — G, 
Ji 

c\o = -9r.g2 + — ( drgi + —dtgi + ^drfi 
91 V 92 fi 

= —drg2 = —F => LJoil = F. 



(70) 



This highlights that the functions F and G defined earlier 
are in fact components of the spin-connection; indeed this 



is how Lasenby et al. ( 1998 1 originally defined them, but 



using geometric algebra. Hence the components of the four- 
acceleration in the tetrad frame for general radial motion 
are 



^i,^+G{vy + Fv%' 



(71) 



If we now specialise to the case where v = u, the 
four-velocity of our observer, then [u^] = [1,0,0,0], and so 
[a'] — [0,G,0,0]. Thus the proper acceleration of our ob- 
server is a = ^ —ayhi — G, which coincides with our earlier 
identification of G as a radial acceleration. Equivalently, the 
invariant force per unit rest mass / (provided, for example, 
by a rocket engine) required to keep the observer in this state 
of motion has magnitude G in the outwards radial direction. 



5.2 Radially- moving test particle 

Let us now consider a particle in general radial motion such 
that its four-velocity components in the tetrad frame are 



= [cosh'i/j(T), sinh?/)(r), 0, 0], 



(72) 



where tp{T) is the particle's rapidity in that frame and r is 
the particle's proper time. Using the tetrad definitions ([7| 
one can show that these components are related to those in 
the coordinate basis by 



t = /i coshi/;(r), 

r = g2 coshi/)(r) + gi sinh'!/;(r). 



(73) 



^ jk=e x[Okej + r „„ek ) 



(67) 



Substituting the coefficients ( 72 1 into the equations in 
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(711, we find that the components of the particle's four- 
acceleration in the tetrad frame are given by 



a — sinh 



iP{t) ^^(r) + G cosh ip{r) + F sinh iP{t)^ , 

iP{t) [i/'(r) + G cosh V'(r) + F sinh 7/>(r)j . (74) 



a — cosh 



Thus, the particle's proper acceleration a = y/—a^ai, and 
hence the invariant force per unit rest mass / required to 
maintain the particle in this state of motion, is 



/ = i>{T) + Gcoshi/'(-r) + Fsm\iip{T). 



(75) 



It is worth considering also the expression for r in terms 
of the force / applied to the particle. Differentiating equa- 



tion ( 73 1 with respect to the proper time of the particle, one 



obtains 

r = iji sinh?/'(r) -f 32 cosh7/;(T) 
-I- V'(t) [gi cosh ?/>(r) + 32 sinh V'(t)] , 



(76) 



which, on substituting for i/^(r) from (751, gives 

r = g\ sinh?/'(r) -I- 32 cosh7/;(T) 
+ [/ - Gcosh7/'(T) - Fsinh7/>(r)] 

X [gi cosh ■!/'(r) + p2 sinh i/;(r)] . (77) 

5.3 Force required to keep test particle at rest 

We now consider the special case of a particle at rest relative 
to the central point mass. In this case, f = and from (731 
it is found that sinh-(/)(r) = — ((72/gi) cosh?/;(T). Moreover, 
since the magnitude of the four-velocity must be unity, we 
deduce that sinh?/'(r) = ±52 1 sj — 92 ■ Hence, the four- 
velocity ( 72 1 may be written 



[51,-32,0,0] 



(78) 



Using these expressions for cosh-(/;(r) and sinh7/;(r) in 
( 75 1, the force per unit rest mass / required to keep the test 
particle at rest is found to be 



1 



9291 - 9291 
\/9x - 9i 



+ G51 - F92 



(79) 



For general radial motion, djdr — tdt+rdr, but in this case 
r = and from (731 and (781 we find i = j\9\l \J 9-^ - 92 ■ 
We thus obtain 



1 



\f9F^^' 



fi9i{92dt9i - 9i9t92) 

2 2 

9i - 92 



G91 - F92 



(80) 



Hence, the required force is determined by the tetrad com- 
ponents themselves and also some of their derivatives with 
respect to t and r (the latter through the functions F and 
G). 

It is of interest to compare our expression for / in ( 79 1 
with the expression ( |77| l for r in the special case of a particle 
released from rest, for which / = and r = 0. In this case. 



using (78 1, we find that, instantaneously, 



9291 - 9291 



92^ 



G91 - F92 



(81) 



Thus, we see that the force / (per unit rest mass) in (79 1 
required to keep a test particle at rest is related to the in- 
stantaneous value of r for a particle released from rest by 



'-92' 



(82) 



The force expression ( |80[ ) in this general form can be 
applied to any spherically-symmetric system for which the 
required quantities can be computed. Below, we will obtain 
the expressions for / in each of our newly-derived metrics. 
Before we consider each metric separately, however, it is 
worth noting that the force ( |80[ ) becomes singular at any 
point where 91^ = 92'. In the case of our newly derived 
metrics, we see from Table [l] that this corresponds to where 



1 



2m 



WW) 



which is valid for A: = and k 
ture density parameter Q.k{t) - 



r'H'{t) = Q, (83) 
±1. Introducing the curva- 



-k/\R{t)H(i)Y 



where fl(t) is the total density parameter, the above condi- 
tion becomes 



(84) 



Comparing this condition with (34 1, we see that one has 
simply replaced R'{t) by l/\Q( t)H (t)]. Consequently, pro- 
vided m ^ l/[3^/Sfl{t)H^{t)], (84 1 has three real roots, one 
of which lies at negative r and the remaining two lie at 



ri{t) and r2{t) given by (35 1, but with R{t) replaced by 
1 / ^^/ fl{t) (t) . In particular, the force / becomes infinite 
outside the Schwarzschild radius r = 2m and inside the 
'scaled' Hubble radius r = 1/ ^n{t)H^(t). 



5.3.1 Spatially-flat universe 

The functions defining the metric for a point mass embedded 
in a spatially-flat universe (fc = 0) are given in Table [l] To 
calculate /, we must also evaluate gi and g2, which are easily 
shown to be 



91 = 0, 



92 = 



rHHt){qit) + l) 



-r^H^it) 



(85) 



Thus the outward force / required to keep a test particle at 
rest relative to the central point mass is 



-rH^{t) 



(1 - ^ - r-^mit)) 



1/2 



+ 



rH\t){q{t) + l)^- 
(l-^-r2iJ2(t)) 



3/2 



(86) 



Comparing ( 86 1 with the expression ( 65 1 giving the in- 



stantaneous value of r for a particle released from rest in 
this spacetime, we see that they indeed obey the relation- 



ship (821. Assuming m <^ r <C 1/H(t), expanding (86 1 in 



small quantities one finds that the zeroth order term is sim- 
ply m/r^, consistent with Newtonian gravity. We are inter- 
ested in where the cosmological expansion parameter H{t) 
enters into the force expression, and we find that it does so 
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at first order in small quantities. Keeping the leading terms, 
one obtains 



+ qit)H\t)r, 



(87) 



which is consistent with the corresponding equation of mo- 



tion ( 62 1 in the Newtonian limit for a particle moving radi- 



ally under gravity. Finally, for comparison against the equiv- 
alent open and closed universe results, note that keeping all 
terms up to second order in small quantities leads to the 
expression 



^ + ^ + ^+ <im\t)r + hq{t) + - ) H\t)m. 



5.3.2 Open universe 

The functions defining the fc = — 1 metric are listed in Table 
[l] and it can be shown that gi and g2 are given by 



91 



92 



rH\t){q(t) + l)\/l-^ + 74y-^i 



Thus the expression for the required outward force is 



(89) 



\rH\i){q{t) + l)(l-''-^ 

+ 



( 



1 2m 
J- — T -f^ 



3/2 



(90) 



In the limit r/R{t) — > 0, this expression reduces to the force 
(1861) in the spatially-flat case, as might be expected. 



Assuming that m ^ r <^ 1/H{t), expanding (90 1 in 
small quantities and keeping only the leading terms at first 



order, the force is found to be given by equation ( 87 1 . Thus 



in this limit the force in an open universe is indistiguish- 
able from that in a spatially-flat universe. The differences 
between the two become visible only when expanding up to 
second order in small quantities: 



2^ 



-I- q(t)H^{t)r + 2q{t) + ^ H^{t)m - ^ 



2i?2(t)- 



(91) 



The only difference between this expression and the corre- 



sponding spatially-flat result ( 88 1 is the extra negative con- 
tribution proportional to m/R (t). This can be interpreted 
as an 'anti-gravitational' force experienced by the test par- 
ticle due to a virtual 'image' mass dragged out to the curva- 
ture scale of the universe. Since R{t) ^ r this term is neg- 
ligible compared to the m/r^ term, and hence the force is 
not signiflcantly altered relatively to the spatially-flat case. 



5.5.5" Closed universe 

The functions deflning the closed universe {k = 1) metric 
are listed in Table [l] In this case, 



91 



92 = 



RHt)^l-^-^-^^-r^HHt) 
rH\t){q{t) + 



r^H^it) 



(92) 



Thus the expression for the required outward force is 

m I {fl-l)r ^H-2^,\ 



/ 



+ 



[rH'it){q{t) + l) 



1 - 



2m i-^ 



■1 2m 
1 Z -BJ 



r^H^it) 



3/2 



(93) 



In the limit r/R{t) — >■ 0, this expression also reduces to the 



force ( 86 1 in the spatially-flat case. 



Assuming that m <C r ^ l/H{t), expanding (93 1 in 



small quantities and keeping the leading terms up to first 
order, the force is once again found to be given by equation 
(87 1. In this limit the force in a closed universe is thus in- 
distiguishable from that in a spatially-flat or open universe. 
Any differences again only become visible when expanding 
up to second order in small quantities: 



m ^ m 



+ 



2r^ 



+ q{t)H\t)r + 2q{t) + - H\t)m + 



2i?2(t)- 



(94) 



Thus we see that the only difference between this force and 



the corresponding result (911 for an open universe is the 
sign of the extra contribution proportional to m/R^{t). In a 
closed universe, this contribution to the force by the 'image' 
mass, which in this case can be viewed as a real object, is of 
the same sign as that due to the original mass. Once again, 
however, since R{t) ^ r this term is negligible compared 
to the m/r^ term, and hence the force is not significantly 
altered relative to the spatially-fiat case. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a tetrad-based method for 
solving Einstein's field equations in general relativity for 
spherically-symmetric systems. Our method is essentially a 
translation of the method originally presented by |Lasenby| 
et al. ( 1998 1 in the language of geometric algebra. 



We use this technique to derive the metrics describing 
a point mass embedded in an homogeneous and isotropic 
expanding fiuid, for a spatially-ffat, open and closed uni- 
verse, respectively. In the closed universe case, the lack of 
the notion of spatial infinity means that considerable care is 
required to determine the appropriate boundary conditions 
on the solution. We find that in the spatially-fiat case, our 
metric is related by a coordinate transformation to the cor- 



responding metric derived by McVittie (19331. Nonetheless, 
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our use of 'physical' coordinates greatly facilitates the phys- 
ical interpretaion of the solution. For the open and closed 
universes, however, our metrics differ from the correspond- 
ing McVittie metrics, which we believe to be incorrect. 

We derive the geodesic equation for the motion of a test 
particle under gravity in our spatially- fiat metric. For radial 
motion in the Newtonian limit, we find that the gravita- 
tional force on a test particle consists of the standard 
inwards force due to the central mass and a cosmological 
force proportional to r that is directed outwards (inwards) 
when the expansion of the universe is accelerating (deceler- 
ating). For the standard ACDM concordance cosmology, the 
cosmological force thus reverses direction at about z ~ 0.67, 
when the expansion of the universe makes a transition from 
a decelerating to an accelerating phase. It is expected that 
this phenomenon has a significant impact on the growth and 
evolution of structure, particularly around z ~ 1; a detailed 
investigation of this is a subject for further study. Nonethe- 
less, we believe this effect is implicitly included in existing 
numerical simulation codes. In our companion paper NLH2, 
we will show how the cosmological force leads to a maxi- 
mum size for a bound object such as a galaxy or cluster, 
and explore the stability of bound orbits around the central 
mass. 

We also used our tetrad-based approach to derive an 
invariant fully general-relativistic expression, valid for arbi- 
trary spherically-symmetric systems, for the force required 
to keep a test particle at rest relative to the central mass. We 
apply this result to our derived metrics in the spatially-flat, 
open and closed cases. To first order in small quantities the 
force in all three cases is the same, but differences become 
visible when considering the second order terms. Interest- 
ingly, we find that the force in an open universe has an 
additional component that may be interpreted as a gravita- 
tionally repulsive term due to an 'image' mass m dragged out 
to the curvature scale of the universe. A similar term, but 
with the opposite sign, is present in the force in a closed uni- 
verse, and results from an image mass at the antipodal point 
of the universe. In both cases, however, the additional com- 
ponents are neglible at distances from the original central 
mass that are small compared with the curvature scale, and 
so in practical terms the force in a spatially-curved universe 
is not significantly different from that in the spatially-flat 
case. 

It must be pointed out that our assumption that the 
background fluid density p{t) is spatially uniform may be 
questionable, since the point mass breaks homogeneity. The 
correct treatment of the background may require p also to 
depend on the radial coordinate r, which would significantly 
complicate the equations and would probably not yield an- 
alytical solutions. It is interesting that McVittie's k — ±1 
solutions do in fact yield background densities dependent on 
both r and t. Nonetheless, since the point mass in our model 
does not occupy any space in the background, and also be- 
cause we have ultimately been interested in m ^ r ^ R{t), 
assuming p = p{t) is a good approximation. Note that we 
will improve our model by accounting for the finite size of the 
central object in a subsequent paper. We will consider only 
spatially homogenous objects for simplicity, but we leave 
the consideration of more general radial density profiles for 
future research. We also note that our approach may be ex- 
tended to systems with accretion onto the central mass by 



the replacement of m — >■ m{t) in equation (121; a full anal- 
ysis of this will also be presented in a future publication. 
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